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Abstract. We suggest a general method for analyzing aggregate insurance claims that arrive 
according to a very general point process, known in the literature as the order statistic point 
process, which includes as special cases the classical compound Poisson and the Sparre 
Andersen models, among many others. We also allow for the process to govern claim sizes 
via a general dependence structure that relates claim sizes to claim and/or inter-claim times. 
The obtained general results are supplemented with special closed-form illustrative formulas, 
that also exemplify the potential for future research in the area. 

Keywords: Aggregate claims; point process; order statistic; order statistic point process; 
OS-point process; claim frequency; claim severity. 



1 



1 Introduction 



We are interested in building a general and encompassing stochastic model for aggregate 
losses that insurance companies accumulate over certain periods of time. Naturally, the 
model is based on a claim-arrival process, which we seek to be a very general point process 
so that well studied models such as the compound Poisson and Sparre Andersen models 
would be included. We refer to, for example, Li (2008), Ren (2008), Leveille et al. (2010), 
and references therein, for results involving the aggregate (discounted) claim process under 
the Sparre Andersen model. 

Our interest in the problem has also been motivated by other intriguing studies such 
as modeUing hurricane losses, as initiated in the works by Garrido and Lu (2004), and Lu 
and Garrido (2005, 2006). These authors assume that claims arrive according to a non- 
homogeneous Poisson process, which can reflect the seasonality and other periodic features 
of claim arrivals. 

In addition to the general claim arrival process, we also tackle situations when claim sizes 
are dependent on preceding inter-claim times, as in the work by Boudrault et al. (2006). Such 
extensions and generalizations facilitate modelling aggregate claims related to, for example, 
volcano eruptions and earthquakes, when longer time intervals between consecutive volcano 
eruptions or earthquakes often result in more severe damages. We tackle even more general 
situations when claim sizes depend on their arrival times, which allows us to incorporate 
features such as the seasonality of certain claim arrivals (e.g., Garrido and Lu, 2004; Lu and 
Garrido, 2005, 2006). 

We next develop some mathematical formalism. Let Ti, T2, . . . denote claim arrival times, 
and let Xi, X2, . . . be the corresponding claim sizes. Let N{t) denote the number of claims 
that arrive up to and including time t > 0. Hence, N{t) is the largest integer such that 

< Ti < • • • < Tjv(t) < t. The aggregate claim size of all the claims up to and including the 
time t > is therefore the sum 

N(t) 

5(i) = ^X„ 

i=l 

where, by definition, S{t) = when N{t) = 0. Hence, the claim-counting process N{t), 

1 > 0, is a right-continuous stochastic process with unit steps at the time instances Tj and 
such that A'"(0) = 0. Denote the inter-claim times by Vi — — Tj-i, with Tq = by 
definition. 

One of the most basic examples of such a claim-counting process is the homogenous Pois- 



2 



son process under which the inter-claim times are independent and exponentially distributed 
random variables. Assuming that the claim sizes are independent and identically distributed 
(i.i.d.), and independent of the claim-counting process N(t), the aggregate claim size S(t), 
t > 0, is a compound Poisson process (e.g., Klugman et al., 2008, Section 11.1.2). In this 
case we have the formulas 

E[S{t)] 



and 



t 

Var[5(t)] 



E[Xi]Ao (1.1) 
= E[X2]Ao, (1.2) 



t 

where Aq is the Poisson parameter. 

More generally, assuming the Sparre Andersen model, for which the inter-claim times are 
i.i.d. random variables with a common arbitrary distribution and independent of the claim 
sizes, we have that E[S{t)] = E[iV(t)]/i and Var[S(t)] = Var[A^(t)]/x2 + E[N{t)]a^, where /x 
is the mean of Xi, and o"^ is the variance of Xi. 

In this paper we consider a very general claim-counting process, which is called the order 
statistic point process or, for short, the OS-point process. It includes such classes of point 
processes as the non-homogenous Poisson process, the linear birth process, and a number 
of other ones (see, e.g.. Berg and Spizzichino, 2000, and Debrabant, 2008). We concentrate 
our attention in this direction of generality and thus deliberately restrict ourselves to only 
calculating the mean E[S'(t)] and the second moment E[5'^(t)] of the aggregate claim Sit), 
thus keeping mathematical complexities within reasonable limits and facilitating a clearer 
and more transparent exposition of main ideas, results and their proofs. If desired, the herein 
presented ideas and techniques can be extended to more complex functionals of S{t) such as 
its Laplace transform, moment generating function, etc. 

Furthermore, we allow each claim size Xi depend on its arrival time Tj as well as on 
the preceding inter-claim time Vi, or equivalently on Vi and Tj_i, or perhaps just on the 
inter-claim time Vi. Such dependence structures reflect our belief that, in some situations, 
the longer a claim does not arrive, the larger the claim is expected to be. This is a common 
situation with claims resulting from earthquakes, hurricanes, and other catastrophic events. 

The rest of the paper is organized as follows. In Section H] we recall basic facts about 
the order statistic point process, introduce some notation, and derive general formulas for 
the mean E[S'(i(:)] and the second moment E[S'^(t)]. In Section [3] we specialize these general 
formulas to the mixed Poisson process, and then even further specialize them by assum- 
ing the dependence structure between the claim sizes and inter-claim times as suggested by 



Boudreault et al. (2006). These special cases, though obviously of independent interest, 
also illustrate how our general results of Section [2] can be utilized in particular situations. 
Section m develops analogous results for the second moment E[5'^(t)] and thus for the vari- 
ance Var[S'(t)] in the case of the mixed Poisson process and the dependence structure of 
Boudreault et al. (2006). Section [5] concludes the paper. To facilitate an easier reading of 
main results, we have relegated most of the proofs to an appendix. 



2 The OS-point process and general formulas 

Following Crump (1975), we say that N(t), t > 0, is an order statistic point process (i.e., 
OS-point process) if for every t > and n > with the positive probability 

7ri,„ = P[iV(t)=n], (2.1) 

we have that, conditionally on the event N{t) = n, the claim arrival times Ti, . . . ,T„ have 
the same joint distribution as the order statistics of i.i.d. random variables ri, . . . , r„ with a 
common cumulative distribution function (cdf) Ft{x), x > 0, such that Ft(t) = 1. 

Crump (1975) has proved that the cdf Ft{x) is related to the process N(t), t > 0, via the 
equation 

, , E\N(x)] , , 

For example, the non-homogenous Poisson process A^(t), t > 0, with the finite cumulative 
intensity function A(x) = \{y)dy, x > 0, is an OS-point process, where X{y) is the 
intensity function. In this case we have that 

Ft{x) = ^, 0<x<t. (2.3) 

Hence, the homogenous Poisson process with a constant intensity function \(y) = \q > 
is an OS-point process with the cdf 

Ftix) = j, 0<x<t. (2.4) 

When conditioned on N{t) = n, the n claim arrival times Tj can be viewed as the order 
statistics of i.i.d. random variables ri, . . . ,r„ with the common [0, t]-uniform distribution. 
When conditioned on N{t) = n, the n inter-claim times Vi = Ti — Tj_i can be viewed as 
[0, t]-uniform spacings. 

The above example gives rise to the class of point processes with so-called property V, 
which coincides with the class of mixed Poisson processes. Later in this paper we shall use 



the latter processes to illustrate general formulas derived in earlier sections. For further 
examples of OS-point processes, their generalizations, and additional references, we refer to 
Berg and Spizzichino (2000), and Debrabant (2008). 

We note that actuaries have found mixed processes particularly useful for reasons such 
as modelling claims that fall into several categories (see, e.g., Klugman et al., 2008, Section 
6.10). For instance, if we want to differentiate between car insurance policies depending on 
the primary driver's gender, age or driving history, it is appropriate to model data in each of 
these categories using different distributions and then mix them in order to provide a model 
for the entire pool of insurance policies. Due to its attractive properties, the mixed Poisson 
process is a particularly popular choice when modelling aggregate claim amounts (see, e.g., 
Rolski et al., 1999, Section 4.3.3). 

Coming back to our general framework, the following theorem lays down foundations for 
calculating the expected aggregate loss E[S'(t)] under various dependence structures between 
the claim sizes Xj, their arrival times Tj, and the previous arrival times Tj_i. 

Theorem 2.1 Assume that N{t), t > 0, is the OS-point process with an absolutely contin- 
uous cdf Ft and density ft. Let the distribution of each claim size Xi depend on the claim 
arrival process only via the times Tj_i and Ti. Then 



The proofs of this and subsequent results are relegated to an appendix due to their 
technical nature. 

We next derive a general formula for the second moment E[S'^(t)], which together with 
the above formula for the mean E[S'(t)] produces general formulas for the variance and the 
standard deviation of S{t). 

Theorem 2.2 Let the conditions of Theorem \2.1\ be satisfied, and let the distribution of each 
product XiXj depend only on the times Ti_i, Ti, Tj^i, Tj. Then 




(2.5) 



oo 




(2.6) 



n=l 
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where 



A,n = n f'E,[Xl\T, = y]My){l - Ft{y)r-'dy 
Jo 



X Mx)f,{v)F;-'(x){l - F,{y)y'-'dxdv 

and 

Bt,n = n{n-l) [ [ E[XiX2|Ti =y,T2 = z]ft{y)ft{z){l - Ft{z)T-^dydz 
Jo Jo 

X f,{y)fAiii)Uz)[F,[w) - Ffy,))'-^(\ - F,{z)r-'dydwdz 



n-1 



fi=2 



X fMft{w)ft{z)F,{yy-\l-Ft{z)r-^-^dydwdz 

n—2 n , 



(i - 2)\(j - i - 2)\(n - j)\ 

i=2 j=i+2 ^ ' ' ^ ' 



X 



i>y 

E [XiXj\Ti_i = x,Ti = y, Tj_i = w, Tj = z] 

Jo Jo Jo 

X ftix)ftiy)ftiw)ftiz)Fi-'ix){Ftiw) - F,iy)y-'-\l - F,{z)r-^ dxdydwdz. 

Theorems 12.11 and 12.21 allow us to calculate the first and second moments, and thus in turn 
the variance and the standard deviation, of the aggregate loss S(t) based on the conditional 
expectations E[Xj|Tj_i, Tj] and E[XjXj|Tj_i, Tj, Tj_i, T,] as well as on the density ft and the 
cdf Ft. We shall provide specific calculations of such quantities in the following sections, 
where we assume that the claim arrival process is the OS-point process with property V. 



3 Property V and the mean aggregate loss 

Following Feigin (1979), we say that a counting point process N{t), t > 0, satisfies the 
property V if it is an OS-point process and the cdf Ft is given by formula (12.41) . Hence, the 
homogenous Poisson process is a counting process with property V. 

In general, Feigin (1979) has proved that a counting point process has property V if and 
only if it is a mixed Poisson process. With L{X), A > 0, denoting the structure cdf of this 
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process, the probabilities ttj „ - which we from now on denote by 7r^„ - take on the form 

r^e-^*rfL(A). (3.1) 
Jo 

For example, when the structure cdf is degenerate at a point Aq > 0, that is, L{X) = l{Ao < 
A}, then 7r^„, n > 0, are the usual Poisson probabilities with the parameter Aq. 

The next theorem is a variant of Theorem 12.11 in the case of the mixed Poisson process. 
To facilitate a straightforward application of the theorem for analyzing models such as the 
one when each claim size Xi depends only on the preceding inter-claim time Vi (Theorem 
13.21 below), we formulate our next theorem in terms of Xi dependent on the pair (Tj_i, Vi). 
Of course, both pairs (Tj_i, Vi) and (Tj_i,Tj) generate same sigma-algebras and thus convey 
same information about the claim size X^. 

Theorem 3.1 Assume that N{t), t > 0, is the mixed Poisson process with a structure cdf 
L{)^), A > 0. Let each claim size Xi depend only on the time Tj_i of the previous claim and 
the inter-claim time Vi = Ti — Tj„i. Then 

E[S{t)] = Y^7rdn f\[X,\T, = y]^^^^dy 

n=l I ^ 

The next theorem is a special case of Theorem 13.11 when each claim size Xi depends 
only on the preceding inter-claim time Vi. This has been a popular model in the actuarial 
literature (e.g., Boudreault et al., 2006, and references therein). In what follows, we use the 
notation 

n 

Q{n\v) = ^Ai{v) with Ai{v) = 'E[Xi\V = v]. 

i=l 

Theorem 3.2 Assume that N{t), t > 0, is the mixed Poisson process with a structure cdf 
L{)^), A > 0. Let each claim size Xi depend only on the inter-claim time Vi. Then 

E[S{t)]= Xe-^''E[Q{Nxit-v) + l\v)]dvdL{X), (3.3) 

Jo Jo 

where Nx is the homogenous Poisson process with the constant rate A > 0. 

The following corollary is a simplification of Theorem 13.21 when the conditional random 
variables Xi\Vi = v have same means. The corollary is slightly more general than equation 
(3) of Boudreault et al. (2006), which we shall formulate in a moment. 
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Corollary 3.1 Let the assumptions of Theorem be satisfied, and let the conditional 
random variables Xi\Vi = v have same means. Then 

E[S{t)]= Xe'^''E[Xi\Vi = v]{X{t-v) + l)dvdL{X). (3.4) 

Jo Jo 

Corollary 3.2 (Boudreault et al., 2006, eq. (3)) Assume that N{t), t > 0, is a homogenous 
Poisson process with a constant rate A > 0. Let each claim size Xi depend only on the inter- 
claim time Vi, and let the conditional random variables Xi\Vi = v have same means. Then 

E[S{t)]= [ Xe-^''E[Xi\Vi = v]{X{t-v) + l)dv. (3.5) 
Jo 

Alternatively, integrating the right-hand side of equation (13. 5p by parts, we have 

E[S{t)] = Xt(E[Xi\Vi = 0] + e-^" (l - ^) dE[Xi\Vi = v]^, (3.6) 
which appears as identity (3) in Boudreault et al. (2006). 

Boudreault et al. (2006) assume that each claim size Xi depends only on the preceding 
inter-claim time V^, and that the conditional random variables Xi\Vi = v are identically 
distributed with their common cdf given by the formula 

P[X, < y\V, = v] = {l- e-^'") He{y) + e-^^H,{y), (3.7) 

where Hi and Hg are two cdf 's, and /3 > is a parameter. Since larger values of Vi result 
in smaller e~^^\ it is natural to view as the cdf of larger claims and Hg as the cdf of 
smaller claims. Under this dependence model, we have the following corollary with Yi and 
Ys denoting random variables with the cdf's and if^, respectively. 

Corollary 3.3 Assume that N{t), t > 0, is the mixed Poisson process with a structure cdf 
L{X), A > 0. Let each claim size Xi depend only on the inter-claim time Vi,, and let the 
distribution of the conditional random variables Xi\Vi be given by formula (13. 7p . Then 

/»oo /»oo 

E[S{t)] = E[Ye] {Xt-^{t,X,l3))dL{X) + E[Ys] 2t(t, A, /3) rfL(A), (3.8) 
Jo Jo 

where 

2t(t, A, /3) = (/3 - ^ + (/3 + A)At) . 

Under the conditions of Corollary 13. 3[ the following four statements hold: 
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1. We have that, when t — )■ oo, 

2l(t,A,/3) A2 



+ 0(1). (3.9) 



t /3 + A 

A closer look at the remainder term o(l) shows that we can apply Lebesgue's dominated 
convergence theorem and obtain the limit 

,„„ = EM [ (^) AdL(A) + EKI [ (^) Adi(A). (3.10) 

2. If there is only one type of claims, that is, Yi and Ys have same distributions as a 
random variable Y , then for every t > and regardless of the value of /3, we have that 



t 



E[Y] / AdL(A). (3.11) 



Note that the right-hand side of equation (13.111) can be written as E[y]E[A], where A 
is a random variable with the cdf L. 

3. If the structure cdf is degenerate at the point Aq, that is, ^^(A) = l{Ao < A}, then from 
equation fl3.10p we have that 

- >^l>'l (^) ^» - >=W ijTx^) (^-^^l 

4. If the structure cdf is degenerate at the point Aq and there is only one type of claims, 
then from equation (I3.1ip we have that, for every t > 0, 

SM=E[F]Ao. (3.13) 
We can view this result as a ' mixed- Poisson analogue' of equation (11. ip . 



4 Property V and the variability of aggregate losses 

In this section we investigate the variability of the aggregate loss S{t) in the case of the 
counting process with property V, that is, when the process N{t), t > 0, is mixed Poisson. 
We choose the variance of S{t) to measure the variability, and for this, given the results of 
the previous section, we only need to derive formulas for the second moment E[S'^(t)]. 

We start with a reformulation of Theorem 12.21 under the assumption that the counting 
process is mixed Poisson and claim sizes Xi are governed by Tj_i and Vi = Ti — Tj_i, instead 
of Tj_i and Tj. 
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Theorem 4.1 Assume that N{t), t > 0, is the mixed Poisson process with a structure cdf 
A > 0. Let each claim size Xi depend on the claim arrival process only via the time 
Ti-i of the previous claim and the inter-claim time Vi = Ti ~ Ti^i, and let the distribution 
of the product XiXj depend on the claim times Tj_i and Tj_i as well as on the inter-claim 
times Vi and Vj. Then 

oo 

E[52(t)]=5^vr^„(A,^„ + 2Sj„), (4.1) 

n=l 

where 



Ai^ = nj\[Xl\T,=y]^^^^dy 



and 



B^^^ = n{n-l) / E[XiX2|Ti = = t^] ^ '^—^dvdy 

Jo Jo 

(j - 3)\in - j)\ I I I 



" 1 Pt f-t j-t-W 



i=3 



(w — yV ^{t — w — vY ^ , , , 
X ^ ^ '- dvdwdy 



n—l I 



(2-2)!(n-«- 1)! Jo Jo Jo 

y'-'^it-y -u-vY~'~^ , , , 
X dvdudy 

71—2 n. 



+ EE7— - 



n! 



/•f: nt—X pt pt—W 

Jo Jo Jx+u Jo 



X 



X* ^(w — X — m)-' ' '^(t — w — v)"- 
X dvdwdudx. (4.2) 

Theorem 14.11 is a general result that allows us to tackle various dependence structures 
between Xi and the pair (Tj_i, Vi). In particular, in the next theorem we present a formula 
for the second moment E[5'^(t)] assuming that each Xi depends only on Vi. Naturally, in 
this special case we get a shorter and more elegant formula than the one in Theorem 14. 1[ 
We use the notation 

n 

Q{n\v) = J2 with Qi{v) = E[Xf\Vi = v] 

i=l 
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and 

T{n\y,v)= J2 My)^Av) with A,{v) = E[X,\Vi = v]. 

l<i^j<n 

Theorem 4.2 Assume that N{t), t > 0, is the mixed Poisson process with a structure cdf 
L{X), A > 0. Let each claim size Xi depend only on the inter-claim time Vi, and let the 
conditional variables Xi\Vi be independent. Then 

POO ft 

E 



poo ft 

[S^{t)]= / / \e'^''E[Q{Nx{t-v) + l\v)]dvdL{\) 
Jo Jo 

/•oo rt i*t — y 

+ / / X^e-^^y+''^E[T{Nxit-v) + 2\y,v)]dvdydL{X). (4.3) 
Jo Jo Jo 

We next specialize Theorem 14.21 even further. 

Corollary 4.1 Assume that N{t), t > 0, is the mixed Poisson process with a structure cdf 
L{X), A > 0. Let each claim size Xi depend only on the inter-claim time Vi, and let the 
conditional variables Xi\Vi = v be independent and have same first and second moments. 
Then we have that 

/•oo rt 

E[S'^{t)] = / / Xe'^''Qi{v){X{t-v) + l)dvdL{X) 
Jo Jo 

J J X^e-^^y+'"^/\i{y)Ai{v)[{X{t-y-v) + 2Y -2^dvdydL{X). (4.4) 

Corollary 14. 1 1 connects our earlier results concerning £[5^(^)1 with the dependence model 
specified by formula (13.71) . This makes the contents of the next corollary. 

Corollary 4.2 Assume that N{t), t > 0, is the mixed Poisson process with a structure cdf 
L{X), A > 0. Let each claim size Xi depend only on the inter-claim time Vi,, and let the 
distribution of the conditional variables Xi\Vi be given by formula ( \3.7\} . Then with Yg and 
Yg denoting random variables with cdfs Hi and Hg, respectively, we have that 

poo poo 

E[S\t)] = E[r/] / (At - 2l(t, A, /?)) dL{X) + E[F/] / 2l(t, A, (3) dL{X) 
Jo Jo 

/oo 
(^(t, A, 0, 0) - 2® (t, A, 0, P) + m{t, A, P, P))dL{X) 

+ 2E[r,]E[y;] (^(t, a, O, /?) - ^{t, X, /3))ciL(A) 

poo 

+ {E[Yg])' ®(t,A,/3,/3)rfL(A), (4.5) 

^0 
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where 2t(t, A, /3) is defined in Corollary and the ^-quantities are as follows: 
^(f X R - ^^^'(^ - I f I f2 2 A + /3 (1 + t (/3 + A))) 

(4.6) 

and 

SftA.O./!) = -^+t^ + t=— + (4.7) 

tfzi/i 5o^/i formulas \4-(^ ^'^^ ( [^. 7| j m ^/ie case /3 = reducing to 

Q3(t,A,0,0) = t^A^ (4.8) 

We next analyze the result of Corollary 14.21 in detail. First, we recall the asymptotic 
formula (13. 9p . Then we check that, when t — )■ cxd, 

Q3(t,A,/3,/3) 4/3A3 A^ , , 

and 

Of course, equation (14.81) gives 

(4.11) 

A closer look at the two remainder terms o(l) in equations (14. 9 p and (I4.10p shows that we can 
apply Lebesgue's dominated convergence theorem and then, using Corollary 14. 2[ we obtain 
the limit 

+ (EW)' r f7^^+^7^T^V^(A). (4.12) 



It follows that, under the conditions of Corollary 14. 2 [ the following four statements hold: 
1. We have that 

- r l'^™]^ + EK|^)\-^.L(A). (4.13) 
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2. If there is only one type of claims, that is, Yi and Yg are distributed as a random 
variable Y, then 



lim = iE[Y]y ydL{X). (4.14) 







Note that the right-hand side of equation f l3.1ip can be rewritten as (E[y])^E[A^], 
where A is a random variable with the cdf L{X). 

3. Combining statements f l3.1ip and f l4.12p . we have that 

lim Y^^l^ = (E[r])2Var[A]. (4.15) 

Note that when the structure cdf L{\) is degenerate at a point Aq > 0, then Var[A] = 
and so statement ( 14.15^ implies that Var[S'(t)]/t^ — )■ when t — )■ oo. 

4. Statements fl3.1ip and fl4.12p help us to establish an asymptotic formula for the variance 
Var[S'(t)] in the case Var[A] = 0. Namely, we have that 

lim Y^^l^ = (E[r])2E[A]. (4.16) 

t—^OD t 

We can view this result as a ' mixed- Poisson analogue' of equation (II. 2p . 



5 Conclusions 

We have demonstrated that the order statistic point process - which is a remarkably general 
process - provides a tractable model for insurance claim arrivals and enables us to calculate 
various quantities of interest. More complex quantities than the herein tackled first and 
second moments of aggregate claims can be studied using techniques of the present paper, 
which can further be extended and generalized to incorporate even more general dependence 
structures between claim sizes and their (inter-) arrival times. In summary, we believe that 
the herein suggested methodology opens up a fruitful direction for research in Ruin Theory 
and beyond. 
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A Appendix: Proofs 



Proof of Theorem 12. IL Using repeated conditioning, we obtain that 

oo C ^ 

E[S{t)]=Y,^,Jj2^[X,mt)=n] 

n=l i=l 
(ft 

= V7ri,J / E[X^\T^=y,N{t) = n]dF,^t^^{y) 
n=i ^ -^0 

+ E[X,\T,_i=x,Ti = y,N{t) = n]dFi_,^i\t^n{x,y)], (A.l) 

._n J Jo<x<y<t ) 



where 



and 



Fi\t,n{x) = 'P[Ti<v\N{t)=nl i>l, (A.2) 



Fi,,-|i,„(x, y) = P[Ti < X, T, < y\N{t) = n], 1 < 2 < j. (A.3) 

Conditionally on N{t) = n, the random variables Ti,...,T„ can be viewed as the order 
statistics ri.„ < ■ ■ ■ < r„.„ of i.i.d. random variables Ti, . . . , r„ with the common cdf Ft given 
by formula (I22])- Hence, = P[ri:„ < x] and = P[ri;„ < x,Tj.,n < y]- 

Consequently (see, e.g., David and Nagaraja, 2003), the density of -Fn^ „ is equal to 

fi\tAx) = nft{x){l - Ft{x)r-^ (A.4) 

for all X G [0,t] and vanishes for all other x. For 2 < i < n, the density of Fi_i^i\t,n is 

/.-M|t,n(x,y) = ^^_^y '^^_^y ft{x)ft{y)Fl-\x){l - Ft{y))-' (A.5) 

for all (x, y) such that < x < y < t and vanishes for all other {x, y). Plugging expressions 
(1A.4P and (lA.Sp into equation (lA.ip gives (12.51) and concludes the proof of Theorem 12.11 □ 

Proof of Theorem 12.21 Conditioning yields 
nS\t)] = J^vr.J Y.E[X^\N{t) = n] +2 5^7r,J ^ E[X,X,\N{t) = n] . (A.6) 

ra=l i=l ^ n=l l<i<j<™ ^ 

We next calculate the two sums in the curly brackets on the right-hand side of equation 
(]A.6p . Starting with the first sum and proceeding just like in the proof of Theorem 12. II but 
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now with instead of Xj, we obtain 

X^E[x2|iV(t) =n]= fE[X',\T, = y,N{t) = n]dF,\t,„{y) 

i=2 

Y.[Xl\T, = y]f,{y){l-F,{y)Y-'dy 

X Mx)Mv)F:-HxU1 - F,{y)r-'dxdy, (A.7) 

where the right-most equation follows from equations flA.4l) and ( lA.SP . As to the sum in the 
second curly brackets on the right-hand side of equation ( ]A.6p . we decompose it as follows 

n 

J2 E[XiX,\N{t) =n]= E[XiX2\N{t) = n] + ^ E [XiX,-|X(t) = n] 

l<i<j<n j=3 

n— 1 n— 2 n 

+ ^E[X,X,+i|Ar(t) =n]+J2Yl E[X,X,|X(t) = n] 

i=2 i=2 j=i+2 

(A.8) 

and then investigate the resulting four summands separately. We begin with the expectation 
E[XiX2|X(t) =n]= jj E[XiX2|ri = y,T2 = z\dFr,2\t,n{y^ 

n{n-l) [ [ E[XiX2\T, = y,T2 = z] 

xft{y)ft{z){l-Ft{z)r-Hydz, (A.9) 



^0 



where the right-most equality follows from formula (lA.Sp with i = 2. Next, we calculate the 
second sum on the right-hand side of equation (lA.SP and have that 

n 

5^E[XiX,|X(t)=n] 



i=3 



5Z / / / E[XiXj|Ti = y,Tj^i = w,Tj = z]dFij_ij\t^n{y,w, z) 



i=3 



X My)ftiw)ftiz)iFtiw) - Ft{y)y-%1 - Ft{z)r-^dydwdz, (A.IO) 
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where we have used the fact that, for 3 < j < n, the density of -Fi,j_ij|t,„ is equal to (see, 
e.g., David and Nae: 2003) 

= _ 3),(^ _ ^y ftiy)ftiw)Mz){FM - F,iy)y-\1 - F,{z)r-^ (A.ll) 

when 0<y<w<z<t and vanishes for all other {y,w,z). As to the third sum on the 
right-hand side of equation ( lA.SP , we have that 



n-l 



^E[X,X,+i|Ar(t) =n] 



j=2 

n-l 

= E 

i=2 

-E (,_2),(,;_,, !), /„/ i E[A-.Y„|r,_,^,.T.^u,.r.«^.] 

X /,(t;)/,(ti.)/,(j)F,(!,)'-=(l - fi(z))»-'-'%<;uidz, (A.12) 

where the latter follows from the fact that, for 2 < i < n, the density corresponding to 
j is equal to (see, e.g., David and Nagaraja, 2003) 

77 I 

/._i,,+i|,„(l/, w, z) = _ ^ _ ^^^ My)ft{w)ft{z)Ft{yy-\l - F^iz))^-'-' (A.13) 

when 0<y<w<z<t and vanishes for all other {y, w, z). Finally, we compute the fourth 
sum on the right-hand side of equation flA.Sp and have that 

n—2 n 

J2 E ^[X^X,\N{t) = n] 

i=2 j=i+2 

= J^J^ ////^ t^*^^' = x,Ti = y, Tj_i = w, Tj = z] dFi^i^ij_ij\t^n{x, y, w, z) 

1=2 j=i+2 -J -J -J -J 
n—2 n I 



-2j=i+2 (^-2)!(j-^-2)!(n-j)! 



X 



j-t r-z r-w r-y 

/ / / / ^[X,X,j\Ti_^ = x,T, = y,Tj_, = w,Tj = z] 
Jo Jo Jo Jo 

X ftix)ftiy)ftiw)ftiz)Fr\x)iF,iw) - Ftiy)y-'-\l - Ftiz)r-^dxdydwdz, (A.14) 

where the right-most equation holds because, for 2 < i < j < n with j — i > 2, the density 
of Fi_i^.jj_ij|t^„ is equal to (see, e.g., David and Nagaraja, 2003) 

f,-i,^,j-ij\t,n{x,y,w,z) = ^. _ 2)!(j - , L 2)!(n - j-)! 

X Fl-\x){F,{w) - F,{y)y-^-\l - F,{zy--^ (A.15) 
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within the region 0<x<y<w<z<t and vanishes for all other (x, y, w, z). 

With the above formulas, we arrive at equation (12. 6p and finish the proof of Theorem 
O □ 

Proof of Theorem 13.21 By assumption, each claim size Xj depends only on the preceding 
inter-claim time Vi = Ti — Tj_i (by definition. To = 0) and thus the expectation E[Xj|Tj_i = 
x, Vi = v] reduces to E[Xi|Vj = v], which is Aj(f ). Hence, equation (13.21) becomes 



1=2 



/o Jo 



Interchanging the order of integration on the right-hand side of equation (1A.16P and then 
noticing that the resulting inner integral is the complete beta function, we obtain 

= f:<E f Mv)(n^-^^^)dv. (A.17) 

n=l i=l V « / 

Since N{t) is the mixed Poisson process with the structure cdf T(A), using formula (13. ip on 
the right-hand side of equation (lA.17p . we obtain 

(^)(^ 



v)\dvdL{X) 



^0 n=0 ^ i=l > 

POO nt 

= / \e-^''Y.[Q[Nx{t-v) + l\v)]dvdL{\). (A.18) 
This completes the proof of Theorem 13.21 □ 
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Proof of Theorem 13.11 Since Ft{x) = x/t and thus = 1/t for aU x G [0,t\, equation 
(12. 5p becomes 



E[Sit)] = f^^nl n [ E[X,\T, = y f ^J^' dy 

n=l I ^ 



j=2 



In the inner integral on the right-hand side of equation ( lA.lOp . we change the variable of 



integration y into v using the relationship y = x + v. Equation (1A.19P becomes 



nSit)] = n / EfXilTi = y Y fj dy 

n=l [ ^ 



^Y.7T^WU—1S\ I I ^[X^\T,^i = x,T., = x + v]'^-^^^-^^dvdx 



{^-mn-^)\Jo Jo ' " ' ' ' t^' 

(A.20) 

Since N{t) is the mixed Poisson process with the structure cdf L{\), we apply formula (13. ip 
on the right-hand side of equation (IA.20p . rearrange terms, and arrive at equation (13. 2p . 
This concludes the proof of Theorem 13.11 □ 

Proof of Corollary ISHl Since E[Xi\V, = v] = E[Xi\Vi = v] for alH > 1, we have that 
Q(n\ v) = nE[Xi|Vi = v]. Since E[Nx{t - v) + 1] = X{t - f ) + 1, equation ([33]) completes 
the proof of Corollary 13. 1[ □ 

Proof of Corollary 13.21 Equation (13. 5p follows from (13.41) by noting that the homogenous 
Poisson process is a mixed Poisson process with the structure cdf L{z) = 1{A < z}. This 
completes the proof of Corollary 13.21 □ 

Proof of Corollary 13.31 Under assumption (13. 7p . we have that 

E[X,\V, = v] = {l- e-f"^) E[Y,] + e'^^'Efn]. 

Applying this formula on the right-hand side of equation (13. 4p . we have that the inner integral 
there is equal to 

E[Ye] [ {l-e-^'')Xe~^"{X{t-v) + l)dv + E[Ys] f e-^''Xe~^'' {X{t - v) + l)dv. (A.21) 
Jo Jo 

The right-most integral of ( ]A.2ip is equal to 2l(t, A,/3), and the left-most one is equal to 

At — 2l(t, A, (3). This completes the proof of Corollary 13.31 □ 
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Proof of Theorem 14.21 Since the claim sizes Xj depend only on the preceding inter-claim 
times Vi and the conditional variables Xi\Vi are independent, quantities and -Bj^„ in 
Theorem 14.11 reduce to the following ones: 



+ y ^ r r^Q,[v f''^^~''~''^^'^ dvdx (A.22) 



and 



^0 



t pt pt—w 



(w — yV ^(t — w — w)" , , , 
X ^ — ^ dvdwdy 



— y — u — 
X dvdudy 

X ^ dvdwdudx. (A.23) 

We now calculate the sum on the right-hand side of equation (]A.22p . Interchanging the order 
of integration, we obtain that 



This implies that 



Kn = / Q^(^) J dv. (A.24) 



i=l 

7L 



We next calculate the integrals that make up the quantity Bj^^. First, interchanging the 
order of integration and utilizing the definition and properties of the beta function (e.g., 
Abramowitz and Stegun, 1972, formula 6.2.2), we obtain 

= Y,n{n-l)l I A,{y)A,{v r \^ ' dvdy. (A.25) 



j=3 



^0 
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In a similar fashion we show that 

n-l 



i-(z-2)!(n-z-l)! io L L ^^H^^M dvdudy 



i=2 



Y^n{n~l)l I Mu)A,^,iv)^ ^„ ^ (A.26) 



=2 



Jo 



Finally, we calculate the double sum on the right-hand side of equation (1A.23I) . Replacing 
w hj w — X — u, we obtain 

n ^1 „f j-t—x r-t j-t—w 

g^E (.-2)!(.-.:2)!(„-.)! i / / /- ^-("'^'W 



J x+Ji J 



X* ^(ly — a; — u)-^ * ^(t — w — f)" , , , , 
X dvawaudx 

j=2 j=j+2 



^^"2 n , pt—u ft—u—v /•t—u—v—w 





(^-2)!(j 


— i — 


2)!(n 










X 










(^-2)!(j 


— i — 


2)!(n 










X 











JO JO 



dvdwdudx 



i=2 j=i+2 



JO 

a;«-2^i-«-2j^^ — X — M — f — w)""-' 



dxdwdvdu 



X / dwdvdu 

Jo 

J2J2n{n-l) / A,HA,(t;)^ ^^rf^rf^. (A.27) 



Using equations (IA.25p - (IA.27p on the right-hand side of flA.23p . we obtain 



JO 



Y.n{n-l) / AM^A^ r \^ ^ dvdy 

j=3 ^0 Jo 

"Zi^ ft ft-u _ ^ _ ti)"~2 

22n{n-l) I I A,{u)Ai+i{vy ^ dvdu 



i=2 



JO 



n ,t rt-u {t-U-vT-^ 



+ E E ^(^ - 1) / / T^r^d^d^ 

i=2j=i+2 ^0 Jo ^ 

n{n-l)f f A.(y)A,(t;) ^^~^~'^^" \ vdy. (A.28) 

^0 i<i<j<„ ^ 
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We have from dH]), dKMf . and COSj) that 

POO pt 

= / \e-^''E[Q{Nxit-v) + l\v)]dvdL{X) 
Jo Jo 

POO pt pt-y f oo -vn/i _ „, o,)n-2 ^ 

+ / |E (^%)! e-"'^(^ll/.^0|^^^y^^(A). (A.29) 

In the above calculations, we have obtained the expectation E[Q{Nx{t — v) + 1\ v)] inside 
the double integral using analogous considerations as those in equation fl3.3p . Similarly, we 
calculate the sum inside the triple integral on the right-hand side of equation (]A.29p . This 
leads us to equation (14. 3p . The proof of Theorem 14 .21 is finished. □ 

Proof of Corollary 14. li Since Oi(f ) = Oi(f ) and Ai(t>) = Aj(t>) for alH > 1, we have 

E[e{Nx{t -v) + l\v)] =E[Nx{t -v) + l]ei{v) 



and 

E[T{Nx{t -v) + 2\y, v)] =E[{Nx{t - y - v) + 2){Nx{t - y - v) + l)]Ai{y)A,{v) 

= ((A(t -y-v) + 2Y-2^ Ai(y)Ai(t;). 

Using these formulas on the right-hand side of equation fl4.3p . we complete the proof of 
Corollary O □ 

Proof of Corollary 14.21 Using formula (13 .Tp . we obtain from equation (14. 4p that 
E[S\t)] = P j' Ae-^^( (1 - e-P'') E[F/] + e-^''E[Y^-\^ {\{t - v) + l)dvdL{X) 

oo 



+ / «B(t,A)rfL(A), (A.30) 







where we have used the notation 



X ( (1 - e-^^) E[Ye] + e-^^E[n]) ((A(t - y - t;) + 2)' - 2) dvdy. 
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We already know from the proof of Corollary 13.31 that 

\e-^^ ( (1 - e-^") E[r/] + e-^"E[y;2]) (A(t - v) + l)dv 

= E[r/] (At - 2l(t, A, + EK2]2i(t, a, P). (A.31) 

This givens the first two summands on the right-hand side of equation (14. 5p . It remains to 
calculate the right-most integral of ( ]A.30p . For this, we rewrite 23 (t, A) as follows: 



A) =(E[r,])'A' ^ e-^y (l - 6'^^) " e"^" (l - e~^") {{\{t - y - w) + 2)' - 2) dvdy 
+ 2E[F^]E[F,]A' ^ e-^^ (l - e"^^) ^ e-^^e"^" (^(A(t - y - t;) + 2)' - 2) c/t;rf?/ 



+ (E[F,])^A2 e-^J'e-^^ ^ ^ e-^^e"^'' [{\{t - y - t;) + 2)^ - 2) 

To simplify the presentation, we next introduce a general notation encompassing the OS- 
quantities in formulas (14.61) . (14.71) . and (14.81) : 

A, e. A) = A^ ^* e-(^+®)^ ^* " ^-{^+^)- (^(A(t - y - t;) + 2)' - 2) 

Note the symmetry: *B(t, A, 0,A) = A,A,0). After a somewhat tedious checking of 
formulas (14.61) and (14. 7p . we obtain that 



^(t. A) = {pyY^Y[^{t, A, 0, 0) - m{t, A, 0, /3) + A, /3, 

+ 2E[F,]E[n] (® (t. A, 0, /3) - ^(t. A, /3, /3) j + (E[n])^^(t, A, /3, 

This finishes the proof of Corollary 14. 2[ □ 
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